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ABSTRACT 

Using the stellar population synthesis, we model the stellar contribution for a sample of 110 double- 
peaked broad-lines AGNs from the Sloan Digital Sky Survey (SDSS). The stellar velocity dispersions 
(<7*) are obtained for 52 double-peaked AGNs with obvious stellar absorption features, ranging from 
106 to 284 km s" 1 . We also use multi-component profiles to fit [O III]AA4959, 5007 and H/3 emission 
lines. Using the well-established Mbh — <J* relation, the black hole masses are calculated to range from 
1.0 x 10 7 to 5.5 x 10 8 Mq, and the Eddington ratio from about 0.01 to about 1. Comparing with 
the known -Rblr — L relation, we can get the factor /, which indicates BLRs' geometry, inclination 
and kinematics. We find that / far deviates from 0.75, suggesting the non-virial dynamics of broad 
line regions. The peak separation is mildly correlated with the Eddington ratio and SMBH mass 
with almost the same correlation coefficients. It implies that it is difficult to detect obvious double- 
peak AGNs with higher Eddington ratios. Using the monochromatic luminosity at 5100A to trace 
the bolometric luminosity, we find that the external illumination of the accretion disk is needed to 
produce the observed strength of Ha emission line. 

Subject headings: galaxies: active — galaxies: nuclei — black hole physics — accretion, accretion disks 



1. INTRODUCTION 

Double-peaked broad Balmer emission profiles have 
been detected in about 200 active galactic nuclei (AGNs) 
(Eracleous & Halpern 1994, 2003; Strateva et al. 2003, 
2006). A systematic survey of 110 sources (mostly broad- 
line radio loud AGNs at z < 0.4) by Eracleous & Halpern 
(1994, 2003) suggested that about 20% sources show 
double-peaked broad Balmer lines. They found some 
characteristics of double-peaked AGNs: large full width 
at half maximum (FWHM) of Ha line, ranging from 
about 4000 km s _1 up to 40000 km s _1 (e.g. Wang 
et al. 2005); about 50% starlight contribution in op- 
tical continuum around Ha; large equivalent widths of 
low- ionization forbidden lines; large [O I]/[0 III] ra- 
tios. Strateva et al. (2003) found about 3% of the 
z < 0.332 AGNs in the Sloan Digital Sky Survey (SDSS) 
are double-peaked AGNs, the fraction is smaller than 
that in radio-loud AGNs sample. Double-peaked lines 
have also been detected in some low-ionization nuclear 
emission line regions (LINERs), e.g. NGC1097, M81, 
NGC4450, NGC4203, NGC4579 (Storchi-Bergmann et 
al. 1993; Bower et al. 1996; Shields et al. 2000; Ho et al. 
2000; Barth et al. 2001). It remains a matter as debate 
of the origin of the double-peaked profiles. 

There are mainly three models to interpret the origin 
of the double-peaked profiles: the accretion disk (Chen & 
Halpern 1989; Eracleous & Halpern 1994, 2003; Gezari et 
al. 2007); the biconical outflow (Zheng et al. 1991; Aba- 
jas et al. 2006), and an anisotropic illuminated BLRs 
(Goad & Wanders 1996). More recently, observational 
test and physical consideration preferentially suggested 
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that double-peaked profiles originate from the accretion 
disk within a radius from a few hundreds R g to about 
thousands R g (R g = GM^/c 2 , Mbh is the SMBH mass) 
(Eracleous & Halpern 1994, 2003; Eracleous et al. 1997; 
Strateva et al. 2003, 2006; Gezari, et al. 2007). In or- 
der to interpret the sparsity of double-peaked AGNs, the 
origin of the single-peaked lines from accretion disk have 
been discussed (Eracleous & Halpern 2003 and the ref- 
erences therein): larger out radius of the line-emitting 
accretion disk; face-on accretion disk; and the accretion 
disk wind. Very few of the double-peaked high-ionization 
line profiles (e.g. CIV) is due to that these high ioniza- 
tion lines are thought to arise in a wind, not in the disk. 

The masses of central supermassive black holes 
(SMBHs) can provide an important tool to understand 
the physics of double-peaked AGNs if we reliably esti- 
mate them (e.g. Lewis & Eracleous 2006; Lewis 2006). 
During the last decade, there is a striking progress on 
the study of central supermassive black holes. The stel- 
lar and/or gaseous dynamics is used to derive the SMBHs 
masses in nearby inactive galaxies. However, for AGNs, 
this method is very difficult because nuclei outshine their 
hosts. Fortunately, we can use the broad emission lines 
from BLRs (e.g. H/3, Ha, Mg II, CIV) to estimate SMBH 
masses in AGNs by the reverberation mapping method 
and the empirical size-luminosity relation (Kaspi et al. 
2000,2005; Vestergaard 2002; McLure & Jarvis 2002; Wu 
et al. 2004; Greene & Ho 2006a). There is a scaling fac- 
tor with larger uncertainty, which is due to the unknown 
geometry and dynamics of broad line regions, BLRs (e.g. 
Krolik 2001; Collin et al. 2006). Nearby galaxies and 
AGNs seem to follow the same tight Mbh — relation, 
where <r* is the bulge velocity dispersion at eighth of 
the effective radius of the galaxy (Nelson et al. 2001; 
Tremaine et al. 2002; Greene & Ho 2006a, 2006b), al- 
though it remains controversial for narrow-line Seyfert 
1 galaxies (NLSls) (Mathur et al. 2001; Bian & Zhao 
2004; Grupe & Mathur 2004; Watson et al. 2007). Lewis 
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& Eracleous (2006) derived the black hole masses from 
a* through fitting absorption lines of the Ca 11 triplet 
(AA8498, 8542, 8662) for 5 double-peaked AGNs. There- 
fore the determination of the black hole mass from in- 
dependent method of reverberation mapping is an useful 
prober to explore mysteries of the double peaked AGNs. 

Significant stellar contribution in double-peaked AGNs 
makes the measurement of c* possible and reliable. Here 
we present our results on the <7» measurements for the 
sample of 110 double-peaked AGNs from the Sloan Digi- 
tal Sky Survey (SDSS) (Strateva et al. 2003). In section 
2, we introduce our working sample selected from Strat- 
eva et al. (2003). Section 3 is data analysis. We present 
the calculations of the SMBH mass and the Eddington 
ratio, and discuss their errors in Section 4. Section 5 is 
contributed to the discuss of the BLRs in double-peaked 
AGNs. Section 6 is the relation between the peak sepa- 
ration and Mass/Eddington ratio. Section 7 is the En- 
ergy budget for double-peaked AGNs. Our conclusion is 
given in Section 8. The last section is our conclusion. 
All of the cosmological calculations in this paper assume 
H = 71 km s" 1 Mpc -1 , tt M = 0.27, and Q A = 0.73. 

2. SAMPLE 

Strateva et al. (2003) presented a sample of double- 
peaked AGNs (z < 0.332) selected from SDSS. They used 
two steps to select candidates: (1) they selected the un- 
usual ones from the symmetric lines using the spectral 
principal component analysis (PCA) (Hao et al. 2003); 
and (2) they fitted the Ha region with a combination of 
several Gaussians, and only selected AGNs better fitted 
by two Gaussians. From the profiles of the broad compo- 
nents, several fitting parameters including the positions 
of the red and blue peaks (A rc d, Abiue), the corresponding 
peak heights (i? re d, #biue), FWHMs were given in their 
Table 3. They investigated the multi- wavelength prop- 
erties of these double-peaked AGNs and suggested that 
Eddington ratios could be large in these SDSS double- 
peaked AGNs (Strateva et al. 2006). 

SDSS spectra cover the wavelength range of 3800-9200 
A with a spectral resolution of 1800 < R < 2100. The 
fiber diameter in the SDSS spectroscopic survey is 3" on 
the sky. At the mean redshift of 0.24 in the sample of 
Strateva et al. (2003), the projected fiber aperture diam- 
eter is 13.2 kpc, typically containing about 80% of the 
total host galaxy light (Kauffmann & Heckman 2005). 
The stellar absorption features in these SDSS spectra 
provide us the possibility to measure the stellar veloc- 
ity dispersion. We did not apply aperture corrections to 
the stellar velocity dispersions because this effect can be 
omitted for z < 0.3 (Bernardi et al. 2003; Bian et al. 
2006). 

3. DATA ANALYSIS 

There are a number of objective and accurate methods 
to measure cr*, including two main different techniques, 
the " Fourier-fitting" method (Sargent et al. 1977; Tonry 
& Davis 1979), and the "direct-fitting" method (Rix & 
White 1992; Greene & Ho 2006b and reference therein). 
With the development of computing, the " direct-fitting" 
method become much more popular. For SDSS spectra 
with significant stellar absorption features (such as Ca 
H+K AA 3969, 3934, Mg IbAA 5167, 5173, 5184 triplet, 
and Ca IIAA 8498, 8542, 8662 triplet, etc.), <r» can be 



measured by matching the observed spectra with the 
combination of different stellar template spectra broad- 
ened by a Gaussian kernel (e.g. Kauffmann et al. 2003; 
Cid Fernandes et al. 2004; Vanden Berk et al 2006; 
Greene & Ho 2006b; Zhou et al. 2006; Bian et al. 
2006). The SMBH mass can then be estimated from 
the Mbh — <y* relation if cr» can be accurately measured 
from the spectrum of AGN host galaxy. 

Here we briefly introduce the method to measure cr*. 
Adopting the stellar library from Bruzual & Chariot 
(2003), we used the stellar population synthesis code, 
STARLIGHT, (Cid Fernandes et al. 2004, 2005; Bian et 
al. 2006) to model the observed spectrum 0\. It mod- 
els the spectrum M\ by the linear combination of simple 
stellar populations (SSP). The model spectrum M\ is: 



M x (x,M Xo ,A v ,v^a*) = M\ a 



': /,; .A'\\ 



)G(V*,<7*) 

(1) 

where 6 jiA = Lf p (tj, (tj, Zj) is the j th tern- 

plate normalized at Ao, x is the population vector, M\ 
is the synthetic flux at the normalization wavelength, 
r\ is the reddening term, Ay is the F-band extinction, 
and G(u*, er») is the line-of-sight stellar velocity distribu- 
tion, modelled as a Gaussian centered at velocity u* and 
broadened by ct*. The best fit is reached by minimizing 
X 2 , 



X 2 (x,M Xo , Av,v*,a*) 



A=l 



M x ) w x ] 



(2) 



where the weighted spectrum w x is defined as the in- 
verse of the noise in observed spectra. We can obtained 
some parameters including ct* by modelling the observed 
spectrum, which are important for our understanding the 
stellar populations in AGNs host galaxies. 

We download 110 double-peaked spectra from SDSS 
DR5 4 in the sample of Strateva et al. (2003), as well as 
the extinction values at g band. We perform the Galac- 
tic extinction correction to observed spectra by using 
the extinction law of Cardelli, Clayton & Mathis (1989). 
Then, the rest-frame spectra with errors and masks are 
obtained by use of the redshift given in their SDSS FITS 
headers. We normalize the spectra at 4020A and the 
signal-to-noise ratio (S/N) is measured between 4010 and 
4060 A. An additional power-law component is used to 
represent the AGNs continuum emission. We check visu- 
ally our modelled spectra sorted by the equivalent width 
(EW) of Call K A3934 line. It is 'found that the fitting 
goodness depends on the S/N (> 5), the absorption lines 
equivalent widths (EW (Call K) > 1.5 A), and the con- 
tribution of stellar lights (see also Zhou et al. 2006). At 
last, 54 double-peaked AGNs are selected, which are well 
fitted to derive reliable er* from their significant stellar 
absorption. Using the two sample Kolmogorov-Smirnov 
(K-S) test, kolmov task in IRAF, the distributions of the 
luminosity at 5100A in the total sample of Strateva et 
al. (2003) and our sub-sample are drawn from the same 
parent population with the probability of 53%. There- 
fore, this sub-sample is representative of the total sample 

4 http://cas.sdss.org/dr5/en/tools/search/ 
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of Strateva et al. (2003) with respect to the luminos- 
ity at 5100A. Fig. 1 shows a fitting example for SDSS 
J082113. 71+350305.02. The final results are presented in 
Table 1, arranging in order of increasing right ascension. 

Using Interactive Data Language (IDL), we also care- 
fully model the profiles of [O III] and H/3 lines. Since the 
double-peaked profile of the H/3 line and the asymmetric 
profiles of [O III]AA4959, 5007 lines, seven-gaussian pro- 
files are used to model these lines carefully, two broad and 
one narrow components for H/3 plus two sets of one broad 
and one narrow components for [O IIIJAA4959, 5007. 
We take the same linewidth for each component of 
O III AA4959,5007, fix the flux ratio of [O IIL1A4959 to 
O III A5007 to be 1:3, and set their wavelength separa- 
tion to the laboratory value. And we also set the wave- 
length separation between the narrow component of H/3 
and the narrow [O III] A5007 to the laboratory value. We 
do the lines fitting for these 54 stellar-light subtracted 
spectra. For spectra without obvious stellar features, we 
do the emission lines fitting for the extinction-corrected 
rest-frame SDSS spectra without stellar-light subtrac- 
tion. At last we obtain the gaseous velocity dispersion, 
<7 S , from the core [O IIIJA5007 line, the total flux of 
[O III]A5007 line and the monochromatic flux at 5100A. 
Our emission-line profile fitting for SDSS J022014.57- 
072859.30 is shown in Fig. 2. 

4. THE MASSES AND THE EDDINGTON RATIOS 
4.1. The stellar velocity dispersions 

In the synthesis, we focus on the strongest stellar ab- 
sorption features, such as Call K, G-band, and Ca IIAA 
8498, 8542, 8662 triplet, which are less affected by emis- 
sion lines. We put twice more weight for these features 
during the stellar population synthesis. After correcting 
the template and the SDSS instrumental resolution, we 
obtain the value of <7* through the direct-fitting method. 

Cid Fernandes et al. (2005) applied the same synthesis 
method to a larger sample of 50,362 normal galaxies from 
the SDSS Data Release 2 (DR2). Their er* is consistent 
very well with that of the MPA/JHU group (Kauffmann 
et al. 2003), the median of the difference is just 9 km s _1 
. They also gave the <r* uncertainty based on the S/N 
(see Table 1 in Cid Fernandes et al. 2005). Typically, 
the <7» uncertainty based on S/N at 4020A is about: 24 
km s- 1 at S/N=5; 12 km s" 1 at S/N=10; 8 km s" 1 at 
S/N=15. In order to use the typical errors suggested by 
Cid Fernandes et al. (2005), we calculate the S/N and 
the starlight fraction at 4020 A. The S/N is the mean flux 
divided by the root mean square (RMS) of the flux in the 
range 4010 to 4060 A. We also performed the method 
of Cid Fernandes et al. (2005) to compute the S/N by 
using the SDSS error spectrum, the results are almost 
the same. Considering the contribution from featureless 
continuum (FC, represented by a power law), we use an 
effective S/N to show the typical error of <r*, where the 
effective S/N is roughly the S/N multiplied by the stel- 
lar fraction. We use the featureless continuum fraction 
as the up limit of nuclei contribution, because the fea- 
tureless continuum can be attributed either by a young 
dusty starburst, by an AGN, or by these two combina- 
tion (Cid Fernandes et al. 2004). Therefore, the effective 
S/N is the low limit. In our sample, the effective mean 
spectral S/N at 4020A for these objects are 7.2 (see Table 



1). Thus the typical uncertainty in er„ should be around 



20 km 



For 13 objects with effective S/N less than 



5, their er» in Table 1 are preceded by colons. 

Here, we also apply this synthesis method to a sample 
of the local AGNs presented by Greene & Ho (2006a). 
Greene & Ho (2006a; 2006b) performed a research on 
the systematic bias of c* derived from the regions around 
Call triplet, Mglb triplet, and Call H+K, respectively 
(Barth et al. 2002). They argued that the Call triplet 
provide the most reliable measurements of cr» and there is 
a systematic offset between c* from Call K line and that 
from other spectral regions. We use our synthesis method 
to their sample in two manners: one is using whole spec- 
trum, the other is just using partial spectrum between 
3200A and 7500A at the rest frame. We put twice more 
weight for the strongest absorption features of Ca H+K 
AA 3969, 3934, G-band, and Ca IIAA 8498, 8542, 8662 
triplet. We find that the values of c* in these two man- 
ners are similar by performing our synthesis method. By 
using the least-square regression, the best fit between the 
cr* from these two manners (cr^ holc and er* artlal , respec- 
tively) is: o-f hol ° = (5.49 ± 3.29) + (0.95 ± 0.03)of artial . 
The spearman coefficient R is 0.97, with a probability 
of Pnuii < 10~ 4 for rejecting the null hypothesis of no 
correlation. 

Because the SDSS spectral coverage is from 3800A to 
9200A, most of double-peaked AGNs with redshift larger 
than 0.083 will not cover the range of Call triplet. For 
those sources with redshift larger than 0.083, we used the 
above formula to obtain the corrected velocity dispersion, 
i.e. al = (5.49 ±3.29) + (0.95 ±0.03) x er*. The corrected 
velocity dispersion is listed in Col. (8) in Table 1. And 
we find that the corrected velocity dispersion is almost 
the same to the uncorrected one. The largest difference 
is about 10 km s _1 , which is less than the typical error 
of 20km s" 1 . 

Based on the SDSS instrumental resolution, we take 60 
km s _1 as a lower limit of cr* (e.g. Bernardi 2003). Only 
for two objects, SDSS J133433.24 -013825.41 and SDSS 
J214555.03 +121034.17, the measured cr* are below/near 
the lower limit (51 km s _1 and 62 km s" 1 , respectively). 
These two objects are excluded from further analysis. 

The a* value from the partial SDSS spectra is also used 
to compare our result with Greene & Ho (2006a), who fit- 
ted the region around Call triplet directly. In Fig. 3, we 
compared our cr* with theirs. We found that the agree- 
ment is quite good, the cr* difference (cr^ H — £r£ artlal ) dis- 
tribution is -2.1 km s _1 with a standard deviation (SD) 
of 22.7 km s _1 . By using the least-square regression con- 
sidering the errors of cr^ H , the best fit between cr£* H and 

crr rtial is: cr™ = (-16.66 ± 3.56) + (1.09 ± 0.04)crr rtial 
(solid line in Fig. 3). The spearman coefficient R is 0.86, 
the standard deviation is 2.11, with p nu a < 10~ 4 . 

4.2. The Results 

For these double-peaked AGNs with reliable crj, we 
use the Mbh — 0% relation to derive the SMBH mass (e.g. 
Tremaine et al. 2002 and reference therein), 

/ c \ 4.02 

= 10813 (^OOknT^) M - < 3 > 

We calculate the Eddington ratio, e.g., the ratio of the 
bolometric luminosity (ibol) to the Eddington luminos- 
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Fig. 1. — Sample fit of the synthetic population model for SDSS J081700. 40+34556. 34. The top panel is the logarithm of the observed 
spectrum (top black curve) and the synthetic spectra (bottom red curve) (shifted down for clarity). In this top panel, we also show the 
region around Ca H+K AA 3969, 3934 and G-band (left); the region around Mg IbAA 5167, 5173, 5184 triplet, Fe 5270A (middle); the 
region around Ca IIAA 8498, 8542, 8662 triplet (right). The residual spectrum is shown in the bottom panel. 
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Fig. 2.— Sample multi-component fitting of the [O III]AA4959, 5007 lines for SDSS J022014.57-072859.30: composite profile(blue curves); 
narrow components (green curves); broad components (red curves); the residue (bottom panel). 
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Fig. 3. — Comparison of our cr t from the partial SDSS spectra 
with that for a sample of the local AGNs presented by Greene 
& Ho (2006a). The solid line denotes our best linear fit, the red 
dash lines are y = x ± 20 km s _1 . Considering the error bar of 
it* presented by Greeene & Ho (2006a), most of local AGNs are 
located between two red dash lines. 
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Fig. 5. — The BLRs sizes versus the monochromatic luminosity 
at 5100A (where /=0.75 for random BLRs orbits). Objects with 
effective S/N in 4020 A less then 5 is denoted as open stars. The 
black dash line is the empirical size-luminosity relation, Rblr = 
22.3 X [AL A (5100A)/10 44 ergs" 1 ] 69 It - days, found by Kaspi 
ct al. (2005). The black dot line is the empirical size-luminosity 
relation,R B LR = 40 X (AL A (5100A)/10 44 ergs i) - 52 It — days, 
found by Bentz et al. (2006). The red solid line is our best fit with 
the fixed slope of 0.69 (i?=0.29), lower by -0.63 dex respect to dash 
line. 
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Fig. 4. — The histograms of the black hole masses (left) and the 
Eddington ratios (right) for these 52 double-peaked AGNs with 
reliable a%. 



ity (XEdd), where L E dd = 1.26 x 10 38 (M bh /M Q ) erg s -1 . 
We use the monochromatic luminosity at 5100A (XL\ 
at 5100A) to estimate the bolometric luminosity, L bo i = 
c B AL A (5100A), where c B = 9 (Kaspi et al. 2000). These 
results are presented in Table 1. 

For the typical uncertainty of 20 km s _1 for er» = 
200 km s _1 , the error of log rr* would be about 0.05 
dex, corresponding to 0.2 dex for logAfbh- The error of 
log Mbh is about 0.4 considering the error of 0.3 dex form 
the Mbh — c* relation (Tremaine et al. 2002). Richards 
et al. (2006) suggested a bolometric correction factor of 
10.3±2.1 at 5100 A. Therefore, the final Eddington ratio, 
iboi/iEdd) has a large uncertainty, about 0.5 dex. 

In Fig. 4, we show the histograms of the black hole 



masses and the accretion ratios for these 52 double- 
peaked AGNs. The black hole masses range from 1.0 x 
10 7 to 5.5 x 10 s M Q with a mean value of log Mbh /Mq 
7.7(1 ± 0.37. Using the Ha FWHM and the 5100A 
monochromatic luminosity, Wu & Liu (2004) estimated 
the SMBH masses and the Eddington ratios for an assem- 
bled double-peaked AGNs sample. Lewis & Eracleous 
(2006) noted that the BH masses from the Ha FWHM 
are not completely consistent with those from the stellar 
velocity dispersion. Our BH masses derived from a% are 
indeed smaller than those from the Ha FWHM (Wu & 
Liu 2004) by about an order of magnitude. This will lead 
to our Lboi/^Edd larger than that from the Ha FWHM 
by almost an order of magnitude. 

The Eddington ratio has a distribution with mean and 
the standard deviation of log(Lboi/-i'Edd) is — 1.13±0.38. 
It is suggested that the accretion disk in double-peaked 
AGNs is in the mode of Advection Dominated Accre- 
tion Flow (ADAF) (e.g. Eracleous & Halpern 2003). 
When the Eddington ratio is below than the critical 
one i B oi/^Edd ~ 0.00280^ (Mahadevan 1997), the 
ADAF appears, where ao.i = a/0.1 is viscous coefficient. 
For our SDSS sample, all objects have Eddington ratios 
larger than this critical value of 0.0028. The present re- 
sults from Fig. 4 clearly show that these double peaked 
AGNs have accretion disks in the standard regime. 

The black hole masses can be independently tested by 
the relation between the black holes and bulges. The 
Appendix gives details of the test. We find the black hole 
masses are consistent from Mbh — Mb u i ge and Mbh — cr* 
relations. 

5. BLRS IN DOUBLE-PEAKED AGNS 
5.1. The size of BLRs 
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We also calculate the BLRs sizes for these 52 double- 
peaked AGNs using the SMBHs masses derived from crj 
and the Ha FWHM. We firstly transform the Ha FWHM 
to the H/3 FWHM by (Greene & Ho 2005b): 

/ FWHM \ (l-03±0.03) 

FWHMh/3 = (1.07±0.07)xl0 3 ^ Q3 kms 

(4) 

From the SMBH masses derived from the velocity dis- 
persions, we can calculate the BLRs sizes: 



that using mass estimation methods based on the size- 
luminosity relationship and the calibration factor of 0.75 
can lead to an order of magnitude over-estimate of the 
SMBH mass (Wu & Liu 2004). 
If BLRs are disk-like with an inclination of 8, the re- 
flation between the H/3 FWHM and the Keplerian disk 
plane velocity, v p , is given by (Wills & Browne 1986) 



R 



BLR 



M bh «) x 5-123 
/ x FWHM| 



It — days, (5) 



H/3 



where / is the scaling factor related to the kinematics 

and geometry of the BLRs, defined by Mbh = f—^r- 2 -, 
V p is the Keplerian velocity in disk plane. For random 
orietation of BLR cloud Keplerian orbits, / is 0.75. 

In Fig. 5, we plot the BLRs sizes versus the monochro- 
matic luminosity at 5100A (assuming /=0.75 for random 
BLRs orbits). Almost all objects are located below the 
empirical size-luminosity relation (black dash line) found 
by Kaspi et al. (2005). The correlation between the 
BLRs sizes and the monochromatic luminosity at 5100A 
is not too strong. The best fit is shown as red solid line 

in Fig. 5 as i? BL R = 5.1 [AL A (5100l)/10 44 ergs" 1 ] ' 69 
lt-days (R = 0.34, P nun = 0.20), lower by -0.64 dex re- 
specting to the empirical relation found by Kaspi et al. 
(2005). When excluding objects with effective S/N in 
4020 A less then 5 (open stars in Fig. 5), for fixed slope 
of 0.69, the best fit gives almost the same line but with 
a smaller R of 0.28. 

If we used /=0.52 (Table 2 in Collin et al. 2006), logi? 
will increased by 0.16 dex, the BLRs sizes of double- 
peaked AGNs are still deviated from the empirical re- 
lation found by Kaspi et al. (2005) (about 0.48 dex). 
FWHMs of the double-peaked AGNs are about twice as 
broad as other AGN of similar luminosity (Eracleous & 
Halpern 1994, 2003; Strateva et al. 2003). If the double- 
peaked AGNs are not systematically more massive than 
other AGNs, equation 5 suggested that BLRs radii can 
be smaller than for a similarly massive "normal" AGNs 
by a factor of 4. 

We have to point out that here we are not really testing 
the -Rblr — L relation in double peaked AGNs, but the 
comparison with i?BLR — L relation allows us to derive 
the factor /. 

5.2. The factor f, BLR inclinations and Non-virial 

BLRs 

If we use the empirical -Rblr — L relation (Kaspi et 
al. 2005) to derive the -Rblr and the SMBH masses, 
A'-fbh(o^), to do the calibration of the factor / (Onken et 
al. 2004), we find that the distribution of / is 0.179 with 
a standard deviation of 0.171. It is not consistent with 
the value of (/) = 5.5/4 presented by Onken et al. (2004) 
and it is about 1/3 of 0.52 ± 0.13 suggested by Collin 
et al. (2006). Our results suggest that the empirical 
relation between the BLRs sizes and the luminosity at 
5100A does not hold for double-peaked AGNs (see also 
Zhang et al. 2007), otherwise the calibration factor / 
should be as low as 0.179. An important consequence 
of the breakdown in the size-luminosity relationship is 



FWHMh/3 = 2 (V- 



+ U p 2 sin 2 9) 1/2 . 



(6) 



where V r is the random isotropic component. We may de- 
rive the scaling factor as / = l/[4(Vr/U p ) 2 + 4sin 2 #]. Ig- 
noring U r , / = l/(2sin0) 2 , and the minimum of / is 0.25. 
For only ten object with / > 0.25, we can derive the incli- 
nation of 9 by the above formula. The mean inclination is 
56 degrees. It is suggested that double-peaked AGNs are 
not preferentially edge-on (Eracleous & Halpern 1994, 
2003; Strateva et al. 2003), most have an inclination of 
less than 50 degrees. If all these objects were nearly edge- 
on, the obscuring torus would prevent us from seeing the 
broad lines. Collin et al. (2006) suggested that objects 
with large FWHMs, inclination effects are not really im- 
portant. In these 52 double-peaked AGNs, 42 objects 
have smaller / less than 0.25. V r can't be omitted, we 
need to consider the random isotropic component, im- 
plying the non-virial dynamics of BLRs in double-peaked 
AGNs. If we use / to trace the non-virial effect, smaller 
/ means strong non-virial effect. We find / have strong 
correlations with log Mbh and log(iboi/^Edd) (see Fig. 
6). Using the least-square regression, we derive the cor- 
relations are: 

log/ = (-1.63 ± 0.12) - (0.65 ± 0.10) log(L hol /L Edd ), 

(7) 

where R = —0.69, p nu ii < 10 , and 

log/= (-4.31 ±0.97) + (0.44 ± 0.13) logM bh /M Q , (8) 

where R = 0.45 withp nu n = 9.6 x 10~ 4 . When we exclude 
objects with effective S/N in 4020 A less then 5 (open 
stars in Fig. 6), the best fits for / relations between 
SMBH mass and the Eddington ratio give a little larger 
R (—0.75,0.52, respectively). We should note that the 
reason of these strong correlations are that / is derived 
from M bh , Rblr [from AL A (5100l) ], and FWHMh/3- 
These strong correlations suggest that objects with larger 
Eddington ratios have strong non-virial effect. 

3C 390.3 is the only double peaked object with the di- 
rect measurements of the time delay and the host stellar 
velocity dispersion. And it locates in the empirical rela- 
tion between BLRs sizes and the luminosity (e.g. Kaspi 
etal. 2005). By cr* of 240 kms" 1 (Onken et al. 2004), we 
find that / is about 1.3, consistent with the assumption 
of random BLRs orbits in 3C390.3. Its small Eddington 
ratio (Lewis & Eracleous 2006), (2 — 4) x 10~ 2 , is also 
consistent with the virial BLRs dynamics (See Fig. 6). 

6. PEAK SEPARATION, BH MASS AND 
EDDINGTON RATIO 

For double-peaked AGNs, the peak separation, At> = 
Kcd — Abiue, has a large variance ranging from about 1000 
km s _1 to about 10000 km s -1 (Table 3 in Strateva et al. 
2003). In Fig. 7, we show the relations between logAA 
and log(iboi/^Edd)) log Mbh- Using the least-square re- 
gression, we derive the correlation between log Aw and 
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Fig. 6. — The scaling factor / versus the Eddington ratio (left 
panel) and the SMBH masses (right panel). Objects with effective 
S/N in 4020 A less then 5 is denoted as open stars. The dash lines 
are / = 0.25. The correlation coefficient for the relation between 
the peak separation and the Eddington ratio is R = —0.69 and 
that between the separation and the black hole mass is R = 0.45. 




Fig. 7. — The peak separation versus the Eddington ratio (left 
panel) and the SMBH masses (right panel). Objects with effective 
S/N in 4020 A less then 5 is denoted as open stars. The correlation 
coefficient for the relation between the peak separation and the 
Eddington ratio is R = —0.33, which is almost the same for that 
between the separation and the black hole mass (R = 0.35). 



log(Lboi/£Edd) to be: 

log Aw = (3.43 ± 0.08) - (0.17 ± 0.06) log(L bol /L E dd), 

(9) 

where R = -0.33, p nuU = 0.017 (see Fig. 7). For the 
relation between log A A and logMbh, 

log Aw = (2.17 ± 0.55) + (0.11 ± 0.08) log(M bh /M ), 

(10) 

where R is 0.35 with p nu n = 0.011. When excluding 
objects with effective S/N in 4020 A less then 5 (open 
stars in Fig. 7), the best fits for log Aw relations be- 



tween SMBH mass and the Eddington ratio give smaller 
R (—0.22,0.21, respectively). We also do the multiple 
regression for the dependence of Aw on log(Lbol/-^Edd) 
and logMbh, 

log Aw = (2.54 ±0.64) + (0.13 ±0.09) log(M bh /M(01) 
-(0.10 ± 0.09) log(L b oi/iEdd), 
where the R— Square correlation coefficient is 0.15. 

Wu & Liu (2004) also studied this correlation and 
obtained an apparent stronger correlation (R — 0.84). 
Their derived strong correlation is mainly due to the 
very strong correlation between the separation and the 
Ha FWHM (R = 0.88). When the Ha FWHM is fixed, 
they found that the partial correlation coefficient is only 
0.05. In this paper, we use the stellar velocity dispersion 
to estimate the black hole masses, and we also find a 
mild correlation between Aw and i bo i/^Edd, which im- 
plies that the peak separation would be smaller for AGNs 
with higher Eddington ratios. It provides clues to why 
previous double-peaked AGNs have lower Eddington ra- 
tios. 

There is increasing evidence for a disk geometry of the 
BLR (see a review of Laor 2007). If we assume that 
peak separation is created by the doppler shift of the 
movement of a thin annulus and the annulus radius cor- 
responds the location where the self-gravitation domi- 
antes (Bian & Zhao 2002) , we have the radius i?sG / Rg = 
(their eq. 18 in Laor & Netzer 
where R g = 1.5 x 10 13 m 8 cm, m 8 = M bh /1O 8 M , 
m is the Eddington ratio and a is the viscosity. The 
maximum separation of the double peaks in units of A 
under the edge-on orientation to an observer is given 

1 /2 

" = 84.9 a^ /9 m~ 2 / 9 ml /9 A, (12) 



3586m~ 2/9 m 4 / 9 a 2 / 9 



AA = A 



-4. 



\Rsg, 

where we use the peak separation in term of separation 
velocity, 



Aw = 7.76 x 10 J 



-1/9 



km s 



(13) 



where a y .i — a/0.1 and Ha wavelength An = 6563A. 

Considering the uncertainties of the fitting results for 
the peak separation correlations with SMBH mass and 
the Eddington ratio, the slopes are consistent with the 
simple theoretical expectation. We have to stress that 
equation (11) is for the maximum separation (edge-on 
orientation) when the disk structure is given. More so- 
phisticated model is needed for explanations of the de- 
pendence of peak separations and Eddington ratios. The 
scatters in Fig 7 may be caused by different orientation 
of the BLR in the sample. 

We also find no significant correlations between the 
ratio of the red peak height to the blue height and the 
Eddington ratio, the black hole mass, the peak separa- 
tion. Since the Keplerian velocity is much below the light 
speed, the Doppler boosting effects could be hidden by 
complex situations of the BLR. 

7. ENERGY BUDGET 

Here we discuss the energy budget of the line-emitting 
accretion disk. Based on the standard accretion disk 
model, the disk radiation as a function of radius £ is 
(Chen et al. 1989): 



F(0 = 1 x 10 1 



2 M 24 cr 3 o[i-(6/o 



l/2l 



ergs cm 2 s 



(14) 
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Fig. 8. — The ratio of the Ha line luminosity (Lfja) to the line-emitting power (Wdisk) versus the monochromatic luminosity at 5100A 
(log XLx (5100 A)). The solid line is Lft a = Wdiski the dash line is Lhq = 0.2Wdisk> and the dot line is Lho = O.lWdisk- We display the 
histogram of the distribution of iHa/Wdisk (right). 



where mg = M^/IQ^Mq, M24 is the accretion rate in 
units of 10 24 j s" 1 , and the dimensionless radius £100 is 
£/(100i? s ). Using M = L hol /((c 2 ) where the efficiency C 
is 0.1, the gravitational power output of the line-emitting 
disk annulus between £1 and £2 is (Eracleous & Halpern 
1994): 



WW&,6) = 7.7 x L bol 




where L ho \ is in units of erg s _1 . It is noted that Wdi s k 
is independent of the black hole mass, when we use the 
typical radius in units of R g . We use the luminosity at 
5100A to calculate the bolometric luminosity. From the 
work of Eracleous & Halpern (2003) and Strateva et al. 
(2003), the inner radius is about hundreds of i? g , and 
the outer radius is about thousands of i? g . The outer 
radius of line-emitting accretion disk is about near the 
inner position of torus. We can assume typical inner and 
outer radii of £1 = 450i? g and £1 = 3000i? g to calculate 
the energy output for these double-peaked SDSS AGNs, 
i.e. Wdisk = 10~ 1 876 x Lbol erg s -1 . The distribution of 
log(LH Q /Wdisk) is —0.55 ± 0.06 with the standard devi- 
ation of 0.57 (see also Fig. 8). If assuming as much as 
20% of the power was radiated as Ha line (dashed line in 
Fig. 8), our results show that only 36 out of 105 double- 
peaked AGNs would generate enough power to produce 
observed strength of Ha emission. If we adopt the value 
of 10%, more objects (83 out of 105 objects) showed the 
energy problem. It implied that the majority of double- 
peaked AGNs need external illumination of the disk (e.g. 
an inner iron torus or corona) to produce the observed 
strength of Ha line (Strateva et al. 2006; Cao & Wang 
2006). 

8. CONCLUSIONS 

We use the simple population synthesis to model the 
stellar contributions in double-peaked SDSS AGNs. The 
reliable stellar velocity dispersions are obtained for 52 



medium-luminous double-peaked SDSS AGNs with ob- 
vious stellar features. We find that: 1) The black hole 
mass is from 1.0 x 10 7 A/ Q to 5.5 x 10 8 M Q and the Ed- 
dington ratio is from about 0.01 to about 1; 2) The factor 
/ far deviates from the virialized value 0.75, suggesting 
the non-virial dynamics of BLRs; 3) The peak separation 
is mildly correlated with the Eddington ratio and SMBH 
mass with almost the same correlation coefficients, which 
can be interpreted in the doppler shift of thin annulus 
of BLRs created by gravitational instability; 4) Based 
on the line-emitting accretion disk model, we need ex- 
ternal illumination of the accretion disk to produce the 
observed strength of Ha line. In the future, using dif- 
ferent models, we would fit the double-peaked profiles to 
constrain the nature of double-peaked AGNs. We can 
also use the double-peaked AGNs to constrain the BLRs 
origin ( Nicastro 2000; Laor 2003; Bian & Gu 2007). 
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APPENDIX 

THE RELATION BETWEEN THE HOST MASS AND THE SMBH MASS 

For the sample of 52 double-peaked SDSS AGNs, assuming XL\ cx A 5 , AL A (5100A) is translated to AL A (5530A). 
The host luminosity in V band is the value of \L\(5530A) multiplied by the stellar fraction at 5530A. We assume that 
the host luminosity in V band approximates the bulge luminosity ibuige i n V-band. We then use the following formula 
to calculate the bulge mass: log(Mb u i g e/-^©) = 1-18 log(Lbuigc/£o) — 1.11 (Magorrian et al. 1998). Fig. 9 shows this 
bulge mass versus the BH mass from er*. We fit with fixed slope as 1, intercept is 2.93 ±0.05, correlation coefficient is 
0.53. The null hypothesis is less then 10~ 4 . We find that the My Ulgc versus the log a relation is consistent with Fig. 4 
of Faber et al. 1997. Therefore, the SMBH mass from Lbuigo and the Magorrian relation agrees with that from er*. 
When we exclude objects with effective S/N in 4020 A less then 5 (open stars in Fig. 9), for fixed slope of 1, the best 
fit gives that intercept is 2.95 ± 0.05, correlation coefficient is 0.49. 

It is suggested that the mass from the Mbh — relation could underestimate the SMBH mass for the massive 
ellipticals with SMBH mass larger than 1O 9 M (Fig. 2 in Lauer et al. 2007). For our 52 double-peaked AGNs, their 
SMBH masses are all less than 10 9 M© from Mbh — Mb u i ge relation. The present results are less affected by Lauer et 
al's findings. 
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TABLE 1 

Results for 52 double-peaked AGNs. 
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NOTE. — Col. (1): Object name. Col. (2): Rcdshift. Col. (3): Equivalent width of Ca K and the error in units of A. Col. (4): % 2 . Col. (5) fraction of featureless 
continuum component. Col. (6): effective S/N from the S/N at 4020 A and the fraction of featureless continuum component. Col. (7): Uncorrected stellar velocity 
dispersion in units of km s ^ , and values for objects with effective S/N < 5 are preceded by a colon. Col. (8): Corrected stellar velocity dispersion in units of km s . 
For 4 objects with Ca II triplet, it is not need to do the correction, which arc shown as *. Col. (9): Corrected gas velocity dispersion in units of km s _1 . Col. (10): log of 
the monochromatic luminosity at 5100 A in units of ergs s — ^ . Col. (11): log of SMBH masses from corrected stellar velocity dispersion in units of Mq. Col. (12): The 
Eddington ratios. Col. (13): The peak separations in units of km s -1 . 



